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ABSTRACT 



Aims. We use kinematical and chemical properties of 754 CoRoT stars to characterise the stellar populations of the Milky Way disc 
in three beams close the Galactic plane. 

Methods. From the atmospheric parameters derived in Gazzano et al. (2010) with the MATISSE algorithm, we derived stellar distances 
using isochrones. Combining these data with proper motions, we provide the complete kinematical description of stars in three 
CoRoT fields. Finally, we used kinematical criteria to identify the Galactic populations in our sample and study their characteristics, 
particularly their chemistry. 

Results. Comparing our kinematics with the Besancon Galactic model, we show that, within 3<x, simulated and observed kinematical 
distributions are in good agreement. We study the characteristics of the thin disc, finding a correlation that is significant at a value 
of 2-cr between the V- velocity component and the metallicity for two different radial distance bins (8-9 kpc and 9-10 kpc; but 
not for the most inner bin 7-8 kpc, probably because of the uncertainties in the abundances) which could be interpreted as radial 
migration evidence. We also measured a radial metallicity gradient value of -0.097 ± 0.015 dex kpc~' with giant stars, and -0.053 ± 
0.015 dex kpc~'with dwarfs. Finally, we identified metal-rich stars with peculiar high [a/Fe] values in the directions pointing to the 
inner part of the Galaxy. Applying the same methodology to the planet-hosting stars detected by CoRoT shows that they mainly 
belong to the thin disc population with normal chemical and kinematical properties. 

Key words. Galaxy: stellar content, disk, structure, evolution, kinematics and dynamics 



1. Introduction 

Understanding the history of the Milky Way requires a thorough 
study of the Galactic populations. The structure and chemistry of 
the Galactic thin and thick discs have already been extensively 
studied; however, most of these studies are either limited to the 
close solar neighbourhood (the closest 1 kpc, Allende Prieto 
2010; Holmberg et al. 2009, 2007; Reid et al. 2007; Nordstrom 
et al. 2004), they explore in more detail some specific Galactic 
directions (Gilmore et al. 2002; Kordopatis et al. 201 lb), or they 
are devoted to the study of the Galactic thick disc (Feltzing & 
Bensby 2009; Veltz et al. 2008; Soubiran et al. 2008, 2003; 
Siebert et al. 2003). These different studies have improved our 
knowledge of the different Galactic components in the solar 
vicinity but we are far from completely understanding them over 
the whole Galaxy. 

The advent of Gaia (Mignard 2005) will strongly constrain 
the structure and composition of the Milky Way because kine- 
matical and chemical properties will be measured for millions 
of stars. Before Gaia, several large-scale surveys, such as RAVE 
(Steinmetz et al. 2009) and the Gaia-ESO Survey (Gilmore et al. 
2012), are collecting data about the Galaxy, helping us to pre- 
pare to exploit the Gaia results. On the other hand, the CoRoT 
(Convection Rotation and planetary Transits) mission is collect- 
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ing light curves for several thousand stars close to the Galactic 
plane towards two diametrically opposed directions. Therefore, 
a non-negligible by-product of this mission is the study of the 
Galactic structure in the directions observed by CoRoT. 

To prepare and support this mission, massive spectroscopic 
observations have been performed, resulting in good precision 
radial velocities for 1534 CoRoT targets (Loeillet et al. 2008). 
Furthermore, atmospheric parameters, i.e. effective temperature 
( T e ff), surface gravity (log g), global metallicity ([M/H]), and 
a-elements abundances ([a/Fe]), have been determined with the 
MATISSE algorithm for 1227 CoRoT targets in three of the di- 
rections observed by the satellite (Gazzano et al. 2010). These 
data represent a good opportunity to identify and characterise 
the different stellar populations composing these Corot Fields, 
to explore the Galactic structure and chemistry in these direc- 
tions, and to explore radial metallicity and abundance gradients. 

The previous analysis by Gazzano et al. (2010, hereafter 
Paper I) aimed at demonstrating the ability to perform robust 
automated spectral classification. In the present paper, we ex- 
tend the characterisation of the stellar populations in three of 
the CoRoT fields by combining kinematical and atmospheric pa- 
rameters. In Sect. 2, we recall some properties of the three ob- 
served samples and present the kinematics analysis. In Sect. 3, 
we compare and validate our results with the Besancon Galactic 
model (Robin et al. 2003). In Sect. 4, we separate and discuss the 
various stellar populations identified in our spectroscopic sam- 
ple. The properties of the thin disc are analysed in Sect. 5. We 
also suspect the presence of a peculiar population with high and 
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non-standard [a/Fe] values and rather high metallicities, which 
is presented in Sect. 6. The impacts of our study in terms of 
planet population in these CoRoT fields are discussed in Sect. 7. 
Lastly, we give our conclusions in Sect. 8. 

2. Stellar properties in the targeted Galactic 
directions 

2.1. Galactic directions studied 

We used the samples of CoRo T/Exoplanet targets analysed in 
Paper I. These stars are located in three of the CoRo T/Exoplanet 
fields, namely the Long Run Anticentre 01 (LRaOl), the Long 
Run Centre 01 (LRcOl), and the Short Run Centre 01 (SRcOl). 
These observations contain relatively bright stars (J < 15, 
2MASS filter system, see Cutri et al. 2003) located close to the 
Galactic plane towards the Monoceros (LRaOl) and the Aquila 
(LRcOl &■ SRcOl) constellations. 




LRaOl _ 
LRcOl . 
SRcOl. 




Fig. 1. Properties of stars analysed in this study in each of the 
three pointing directions. Top : Distribution of the J magnitudes 
of the spectroscopic sample normalised to the total number of 
targets ; Bottom : Distribution of the J magnitudes of our spec- 
troscopic targets given in percentage of the overall population 
in the corresponding CoRo T/exoplanet field. The red dotted- 
dashed represents LRaOl stars, the blue solid line stands for 
LRcOl stars, and SRcOl is shown with a long-dashed green line. 

As presented in Paper I, the targets were observed during two 
different campaigns with slightly different observation strate- 
gies. The first campaign took place in January 2005, which is 
before the launch of CoRoT, while the second one was done 



more than three years later, when CoRoT was in operation. In 
both cases, the observations aimed at both performing accurate 
Doppler measurements and checking and improving the spectral 
classification of stars in the CoRoT exoplanet fields. The spec- 
tral classification of all the stars within the range of magnitudes 
observable by CoRoT was mandatory to identify dwarfs and gi- 
ants and to further prepare the targets to be selected for CoRoT 
observations in the exoplanet fields. Once the light curves are 
analysed, it is also useful to help the screening out of transit- 
ing systems stellar in origin. Prior to the launch, the estimate of 
spectral types and luminosity classes has been carried out based 
on photometric data from both dedicated ground-based observa- 
tions and existing catalogues (Deleuil et al. 2009). This allowed 
to build the CoRoT input catalogue, Exo-Dat, which provides for 
all stars within the potential exoplanet fields of the mission, the 
astrometry, magnitudes in different bands in the visible and near- 
IR, the spectral type, and luminosity class as estimated from the 
photometry. For our two spectroscopy campaigns, we used sim- 
ilar criteria to select our targets in a given FLAMES/GIRAFFE 
field, based on this input catalogue. As presented in detail in 
Paper I and Loeillet et al. (2008), the highest priorities were 
given to bright F,G, K, and M-dwarfs and subgiants with a r'- 
magnitude less than 15.0 in order to ensure that a signal-to-noise 
ratio decent enough could be reached. The angular separation be- 
tween stars for the positioning of fibers within the field of view 
of the spectrograph imposes additional constrains that require 
not using too strict selection criteria. As a very last priority, we 
thus allow the selection of any type of star fainter than magni- 
tude 15.0 in r'. 

The properties of our targets in the three directions are 
summarised in Table 1. Figure 1 shows the distributions in J- 
magnitude of the observed stars and how they compare to the 
overall stellar population in each CoRoT field. It is clear that 
our sample represents a very small fraction of the total pop- 
ulation in a given field. This is due to the limited number of 
FLAMES/GIRAFFE fields we could observe in spectroscopy 
during the nights that were allocated to our programmes along 
with the observation strategy. The percentage of targets we ob- 
served depends on the field as the stellar populations differ from 
one region to another. While LRcOl and LRaOl which were se- 
lected as prime fields for the CoRoT exoplanet programme, are 
densely populated and present an homogeneous stellar density, 
the SRcOl whose selection was driven by the CoRoT asteroseis- 
mology programme, is poorly populated and very inhomoge- 
neous, with regions of the field affected by a strong extinction. 
This explains why our spectroscopic targets in the SRcOl are 
more concentrated toward bright stars. In contrast, the percent- 
age of spectroscopic targets in LRcOl and LRaOl is distributed 
better over the magnitude range, the lowest percentage for LRcOl 
being due to a much higher stellar counts in this region, com- 
pared to LRaOl . 

2.2. Atmospheric and chemical stellar properties 

In Paper I, we used the MATISSE algorithm (Recio-Blanco et al. 
2006) to derive stellar atmospheric parameters, i.e. the effec- 
tive temperature (T e fj), the surface gravity (log g), the over- 
all metallicity ([M/H]), and the a-enhancement with respect to 
iron ([a/Fe]) 1 , for 1 227 stars from their FLAMES spectra in the 
HR9B configuration. These parameters are affected by several 
sources of uncertainty. The internal uncertainty is the numerical 
uncertainty only due to observational noise in the spectra and 
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Table 1. Properties of the three Galactic directions studied. 



CoRoT field 


Nlbt.Stars 


Nf.s 


ID 


bf) 


J 


LRaOl 


457 


404 


212.2 


-1.9 


[8.7 ; 13.9] 


LRcOl 


555 


286 


37.5 


-7.5 


[9.7 ; 14.6] 


SRcOl 


215 


64 


36.8 


-1.2 


[9.0; 13.8] 


Tot. Sample 


1227 


754 






[8.7 ; 14.6] 



Notes. The columns contain the CoRoT field names, the number of stars 
with MATISSE atmospheric parameters, good kinematics parameters, 
and correct SNR spectra (Final Sample, see Sect. 2), the mean Galactic 
longitude and latitude, and the J magnitude range. 



400 

300 
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Fig. 2. Distribution of the total uncertainty for the four atmo- 
spheric parameters for all the studied CoRoT stars. The largest 
uncertainties in T e ff and log g distributions correspond to the 
few (10%) metal-poor cool giant stars of the sample, whose 
SNR ranges from ~ 10 to ~ 130. 



the MATISSE method. This was estimated in Paper I to analyse 
using MATISSE a grid of interpolated theoretical spectra with 
various signal-to-noise ratios per pixel (SNR). This uncertainty 
corresponds to the self consistency of the parameter estimation 
procedure and how it is affected by noise. It is different from the 
relative uncertainty for which other sources of uncertainty have 
to be taken into account (for instance, uncertainty in the atmo- 
spheric parameters induced by radial velocity uncertainties, bad 
normalisation, and other purely instrumental issues). In Paper I, 
we used the multiple observations of fifty stars to evaluate this 
source of uncertainty. Finally, by comparing our atmospheric pa- 
rameters with several reference libraries, we can estimate an ex- 
ternal source of uncertainty. 

In the present study, we estimated this external uncertainty. 
To that purpose, we used the same reference stars sample as 
in Paper I. This sample is composed of the 118 stars from the 
Elodie 3.1 library (Prugniel et al. 2007), the 90 stars from the 
S 4 N study (Allende Prieto et al. 2004) and the 39 giant stars 
from the study by Santos et al. (2009). We separated these refer- 
ence stars sample into eight subsamples. Each of the three atmo- 
spheric parameters ranges was divided in two narrower ranges, 
at T eff = 5500 K, log g = 3.5 dex, and [M/H] = -0.36 dex. To 
ensure the reliability of our uncertainty estimate, we took ev- 
ery source of uncertainty into account, and calculated the total 
uncertainties, adding quadratically the internal, relative (both re- 
ported in Paper I), and external uncertainties. 



The distributions of the total uncertainty for the four param- 
eters are presented in Fig. 2. The largest external uncertainties 
are found for the metal-poor giant cool stars (8% of the whole 
sample, 79 stars with a median SNR of 25), for which we esti- 
mated total uncertainties of <x Tc(r 244 K, o~\ og g — 0.382 dex, 
C[M/H] — 0.156 dex, and cr[ a /F e ] - 0.06 dex. For the majority of 
the sample, with a median SNR of 23, the total uncertainty in 
the determination of atmospheric parameters is <x Xdr 125 K, 
o"iog g — 0.26 dex, <T[m/h] - 0.17 dex, and cr [Q ,/ Fe ] - 0.1 dex. 
This approach, even if pessimistic, is the safest one because we 
try not to neglect any source of uncertainty. Our total uncertain- 
ties may be overestimated because we assumed the literature pa- 
rameters to be perfect whereas each measurement has associated 
uncertainties, and significant differences can be found between 
different parameter estimations performed by different authors. 

To further test the consistency between these different 
sources of uncertainty, we added Gaussian noise to the spectra of 
the S 4 N sample and recalculated the atmospheric parameters at 
the signal-to-noise ratios 50, 30, 20, 10, and 5. At every signal- 
to-noise ratio, the S 4 N values of the [a/Fe] are recovered within 
the total uncertainty. No correlation in the uncertainties with the 
observed stellar parameters was noticed. As far as the theoretical 
correlations are concerned, Fig. 2.2 shows the error ellipses for 
a set of synthetic spectra, as defined and computed in Sect. 3, 
taking the total uncertainty into account. The correlations shown 
in this figure confirm the results already shown in Paper I that 
the error ellipses are very small. 

2.3. Stellar distances 

To estimate the distance of the stars in our sample, we calcu- 
lated their absolute magnitude from the atmospheric parame- 
ters derived by MATISSE. For that purpose, we used the YYmix2 
interpolation code to generate a set of Yonsei-Yale isochrones 
(Demarque et al. 2004) to derive the absolute visual magni- 
tude (My) following the procedure proposed by Zwitter et al. 
(2010) and implemented by Kordopatis et al. (201 la,b). We note 
that the Yonsei-Yale isochrones do not include any evolutionary 
phase after the Red Giant Branch (RGB), in particular the Red 
Clump (RC) where some of our stars may be lying. Nevertheless, 
we insist that this lack does not have any influence on the de- 
rived stellar distances of our pipeline. Indeed, given the uncer- 
tainties in T e ff and log g (140 K and 0.27 dex, respectively, see 
Paper I), we are not able to reliably distinguish the RGB from the 
RC stars. Nevertheless, we did the exercice of selecting the stars 
with 4500<T eff <5 000 K and 2.5<log g< 3.0 (less than 100 stars 
of the total sample) and modified their absolute J-magnitude to 
meet the one given by Girardi & Salaris (2001) of Mj=-0.8 mag 
(for solar metallicity stars). We then reprocessed the entire anal- 
ysis, and no radical changes have been noticed in the velocity 
distributions or in the metallicity distributions of Table 9. One 
of the reasons for that is the relatively small number of stars 
concerned, and the pipeline we used roughly obtain the good ab- 
solute magnitudes for half of these stars. In addition, we recall 
that the tests done on synthetic data from the Besancon model 
did not show any particular distance bias for any type of star, 
in particular the RC stars. On the other hand, the interstellar ab- 
sorption cannot be neglected because our pointing directions are 
close to the Galactic plane. To estimate this absorption, we cal- 
culated the expected unreddened colours for our sample, from 
the atmospheric parameters derived by MATISSE and the T e ff- 
2MASS colour calibration of Gonzalez Hernandez & Bonifacio 
(2009). We inverted their Eq. (10), keeping the only physical so- 
lution. We used the (J - K s ) colour for the extinction estimate, 
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Fig. 3. Error ellipses computed as the difference between the 
input and the derived parameters of our pipeline, for the set of 
2,000 spectra of Sect. 3, with SNR~ 10. These figures illustrate 
the extent of the correlation of the uncertainties in our parameter 
derivations. The 2-dimensional histograms have been obtained 
by binning by 16 K in effective temperature, 0.025 dex in logg 
and 0.016 dex in [M/H]. The isocontours are for 2, 5, 15 and 30 
stars 



and Cardelli et al. (1989) calibrations to transform colour ex- 
cess into absorption. The stellar distances (D) were then calcu- 
lated from the absolute and apparent magnitudes in the J band. 
Gonzalez Hernandez & Bonifacio (2009) give a validity domain 
for applying their calibrations as a function of the T e ff, [M/H], 
and colours. We rejected 1 14 stars that do not match their valid- 
ity domain. We propagated the lcr uncertainties at every step of 
the procedure ensuring a reliable estimate of the total uncertainty 
on the stellar distance. The resulting D, My, {V - J)o, (J - K s )q 
and A j are reported in Table 2 (electronic form). 

Figure 4-a) shows the distribution of the stellar distance for 
the whole analysed sample. The vast majority is located within 
2 kpc from the Sun. Note that the distribution is not symmet- 
ric, presenting a tail towards further distances, up to 6 kpc. 



fl 0.25 

° 0.20 

<" 15 

2 o.io 



o 0.05 



100 1000 10000 

D [pc] 





70 


1 b) 

i + + + + 




60 


1 ♦ 




50 


+ + + + 

L _ A ^ + ^^r^^; 4 




40 


b 


30 






20 







100 



1000 
D [pc] 



10000 




2 4 6 8 10 12 14 16 
D (n) [pc] 

Fig. 4. a) Distance distribution of the CoRoT stars, b) Relative 
uncertainty on the distance estimation as a function of the dis- 
tance {bottom). The dotted line represents the 50% threshold, c) 
Comparison of our distance estimate with parallaxes for the S 4 N 
sample (Allende Prieto et al. 2004). Note that the uncertainty on 
the distance determination from parallaxes is about 1% (hence 
not always represented). 



When dealing with distances, it is instructive to examine the 
relative uncertainty distribution. Figure 4-b) presents the rela- 
tive uncertainty on the stellar distance as a function of the dis- 
tance. It shows that the typical uncertainty on the stellar dis- 
tance is ~ 30%. Only ~ 7% of the stars show uncertainty 
greater than 50%, and they have been discarded in what follows. 
Furthermore, we checked our distance determination by comp ar- 
ing it with HIPPARCOS parallaxes (Perryman & ESA 1997) for 
the stars in the S 4 N sample (Fig. 4-b) 

We used the atmospheric parameters given in Paper I and the 
procedure described §2.2 for computing the total uncertainty, 
assuming an internal uncertainty of <x X |I = 50 K, <x log g = 
0.08 dex, CT[ M / H ] = 0.05 dex, and o-[ a /p e ] = 0.02 dex (cor- 
responding to SNR=20, see Paper I). The correlation between 
the two distance determinations is 0.88, and no bias is detected 
(ctt-d = 1-3 pc), which validates our procedure. We therefore 
conclude that our stellar distance estimate is realistic. We also 
note that no tests could be done on datasets of observed distant 
giant stars with known parallaxes, because no spectra of such 
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stars exist in the wavelength domain and at the resolution at 
which we observed. Nevertheless, tests done in synthetic spectra 
(see Sect. 3) have shown that we recover properly the distances 
for this kind of stars, even at a SNR~10. 

2.4. Velocity components 

The kinematical data needed for the current study are the 
barycentric radial velocity (V ra d, measured by Loeillet et al. 
2008; Gazzano et al. 2010) and the proper motions (fiufib) where 
I and b are the Galactic longitude and latitude, respectively. 
The proper motions were extracted from the PPMXL catalogue 
(Roeser et al. 2010). These authors showed that the UCAC3 
(Zacharias et al. 2010) proper motions are less reliable northern 
than 5 = -20°, which is the case for our observations, and their 
uncertainties are generally low (< 07, >^ 4 mas yr~'). The cross- 
match of the initial 1227 stars from Gazzano et al. (2010) with 
the PPMXL catalogue resulted in a sample of 1074 stars, i.e. 
~87% of the initial sample described in Table 1. We removed all 
the stars with high uncertainties on the proper motion (cry, > 10 
mas yr~') and, following Roeser et al. (2010) recommendations, 
we also removed very high proper motion stars (fi a or fig > 150 
mas yr _1 ) and stars having too few observations in the PPMXL 
catalogue (n b s < 3). We also removed from our sample all stars 
with an uncertainty in radial velocity greater than 1 kms -1 (the 
average uncertainty being < o> riil >^ 0.3 kms -1 ). This cri- 
terium on the radial velocity is important not only in order to get 
small uncertainties on the derived heliocentric velocities (see be- 
low), but also to derive reliable atmospheric parameters from the 
stellar spectra, hence good estimations of the distances. Indeed, 
as shown in Gazzano et al. (2010), for errors larger than half a 
pixel on the Doppler correction of the spectra, the derivation of 
the stellar atmospheric parameters degrades. 

Finally, to interpret our results and characterise the stellar 
populations in these Galactic directions, we need to clean up 
our sample by introducing a cut in SNR. For that purpose, we 
did not consider all the spectra with SNR < 10 in the following 
because their atmospheric parameters might be less reliable, 
and because for these amounts of noise the correlation in errors 
might become considerable. This SNR cut removes 25% of the 
total initial sample of 1227 stars. Combined with the previously 
cited quality criteria, the final sample results in 754 stars 
corresponding to ~ 62% of the initial sample, i.e. 404 stars in 
LRaOl, 286 stars in LRcOl, and 64 stars in SRcOl (see Table 1). 



where k — 4.74047 kms , and (jui,fib) are the true proper mo- 
tions, i.e. not projected on the sky. For that purpose, we con- 
verted the proper motions (fi a ,jj.g) of the Roeser et al. (2010) 
catalogue into the Galactic coordinates system by using classi- 
cal relations, (a p = l2 h 49 m ,S p = 27.4°) as the equatorial co- 
ordinates of the Galactic pole and Iq = 123° as the origin of 
longitudes. 

We remind the reader that in our convention the U veloc- 
ity is directed towards the Galactic centre, the V towards the 
Galactic rotation direction, and the W up towards the Galactic 
north pole. We propagated the lcr uncertainties estimated in 
the previous sections on the space velocity components and the 
Galactocentric coordinates. Clearly, the dominant source of un- 
certainty is the stellar distance, compared to the radial veloc- 
ity and the proper motions. All these quantities are presented in 
Table 3. 

3. Comparisons with the Besancon Galactic model 

We used the Besancon Galactic model (hereafter BGM, see 
Robin et al. 2003) to simulate the three pointing directions de- 
scribed in Sect. 2. These simulations allowed us to test our ob- 
servations in the context of the canonical scenario for Galactic 
structure and chemistry. Besides, by comparing our observables 
to the simulation results, we can check that the observational se- 
lection biases are correctly taken into account. 

For each given field, we performed a BGM request with the 
faint limit in J magnitude, the mean Galactic longitude and lat- 
itude, a solid angle of five square degrees, which is compatible 
with the size of each FLAMES pointing direction, and no extinc- 
tion law. Each request provides a sample of simulated stars with 
their intrinsic properties (absolute magnitude, effective tempera- 
ture, gravity, age, metallicity [Fe/H], U,V,W velocities computed 
without uncertainties), and the corresponding observables (ap- 
parent magnitude, colours, proper motion, radial velocity, dis- 
tance to the Sun, and interstellar extinction). 

The extinction being poorly constrained in the Galactic 
plane, we decided to apply a posteriori to the BGM stars an ex- 
tinction law fitted to our data. We chose to model the extinction 
by fitting 

A j - a + d exp(cD - b) (7) 

0.8 1 1 1 j 1 1 1 1 



The combination of stellar distances and coordinates allows 
us to calculate the Cartesian coordinates of our sample stars 
whose origin is at the centre of the Galaxy: 



Xgc 
Yqc 
Zgc 



Xq 

Yq- 
Zq- 



■ X = Xq - D cos b cos / 
Y — Yq - D cos b sin / 
Z = Zq + D sin b 



(1) 
(2) 
(3) 



where Xq ^ 8.5 kpc, F0 ^ kpc, Zq = 15 pc, and (X, Y,Z) are 
the heliocentric coordinates. We also computed the correspond- 
ing space velocity components relative to the local standard of 
rest 

U = Vrad cos b cos I - kDfib cos I sin b - kD/ij sin I cos b (4) 
V = V^d cos b sin I - kDpi, sin I sin b + kDfi\ cos / cos b (5) 
W = Vrad sin b + kD^ cos b (6) 




D (kpc) 

Fig. 5. Derived extinction law (red line) from the calculated ab- 
sorptions and distances in the LRaOl direction. The typical error 
bar is represented for 500 pc, 1 kpc, 2 kpc, and 4 kpc. 
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where D is the stellar distance and Ay is the absorption in the 
/ band (both presented in Table 2). The resulting a, b, c, and 
d coefficients are given in Table 4 for each field. As boundary 
conditions, we forced these laws to be null at zero distance and 
equal to the Schlegel et al. (1998) values of the absorption at 
very large distances. Figure 5 illustrates this fit in the LRaOl 
direction. 



Table 4. Adopted coefficients for the fit of the assumed extinc- 
tion law applied to the BGM data and number of simulated stars 
in the BGM. 





a 


b 


c 


d 


Nb BGM 


LRaOl 
LRcOl 
SRcOl 


0.671 
0.380 
2.430 


0.850 
1.205 
-0.314 


-0.358 
-3.348 
-0.472 


-1.569 
-1.155 
-1.516 


26763 
100 352 
397 730 



To compare our sample to the BGM, we need to take the ob- 
servational selection criteria into account by biasing each BGM 
request on which we applied our extinction law. This consists in 
reproducing the distributions of the infrared colour (J - K s ) and 
magnitude J (shown in Fig. 1). This procedure is statistically ro- 
bust for each of the three fields since the BGM requests contain 
a sufficient enough number of stars (> 10 4 , see Table 5). 

3.1. Self-consistency check of the distance calculation 
method 

We took advantage of the BGM request in the LRcOl direction 
to test our stellar distance estimates, derived from the determina- 
tion of atmospheric parameters, using simulated FLAMES spec- 
tra. To that purpose, we selected 2 000 stars randomly in this 
simulated sample and interpolated the corresponding theoret- 
ical spectra, from the grid used for the learning phase of 
MATISSE (see Paper I), hence with the same observational set- 
ting (wavelength range and resolving power) as our observed 
spectra (see Fig. 6 left panel). We added gaussian noise to these 
synthetic spectra and analysed them with MATISSE in the same 
way as our observations (Paper I). At SNR=10, we obtained the 
HR diagram presented in the right hand part of Fig. 6, which 
is not distorted but only scattered. The agreement between the 
input and recovered diagrams is very good. 

To study the uncertainty on the retrieved stellar parameters, 
introduced by MATISSE, we calculated the bias and standard de- 
viation for several combinations of the stellar parameters (see 
Fig. 7). The bias is systematically very small and lower than the 
dispersion for any type of stars. We found values of the internal 
uncertainty that are compatible with Paper I. The standard de- 
viation of the three parameters depends strongly on the effective 
temperature and the metallicity: cool and metal-rich stars present 
lower <x values than hot and metal-poor stars. 

Applying the entire procedure described in Sect. 2.3, we thus 
conclude that our procedure does not introduce any internal bias 
in the determination of stellar distances, even at low SNR, as 
illustrated in Fig. 8. We also checked that, given the derived 
distances with our pipeline and the BGM proper motions and 
radial velocities, we are able to perfectly recover the velocity 
components (U, V, W) with Eqs. 4-6. This ensures good con- 
sistency between our results and BGM simulations. We insist 
that this test complements Sect. 2.3, since it also shows that our 
pipeline derives accurate distances for the distant giant stars and 
that the mild correlation in the derived atmospheric parameter 
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Fig. 8. a) Derived distances from the MATISSE parameters using 
the interpolated theoretical spectra having a SNR equal to 10 
versus the distance computed by the BGM. b) Quantiles of the 
relative difference between our stellar distances and the BGM 



uncertainties does not affect the final positions or velocities (see 
Sect. 3.3). 

3.2. Stellar distance distributions comparison 

Comparison of our stellar distance distributions with those pre- 
dicted by the BGM is an important verification to understand 
which regions of the Galaxy were observed and how our selec- 
tion bias limits the interpretation of these data. We, therefore, 
tried to estimate the effects of the parameter uncertainties on the 
distance distributions. For that purpose, we considered, for ev- 
ery star of the sample, that its distance can be represented by a 
Gaussian distribution centred on the derived distance, whose cr 
is the uncertainty on its stellar distance {cr D ). We then performed 
2 000 Monte-Carlo realisations by randomly picking distances in 
this distribution and we estimated a mean distribution, and an er- 
ror bar associated to every bin. To apply the same procedure on 
the BGM simulated stellar distances, which are provided with- 
out associated uncertainties, we added, for all the stars farther 
than 250 pc, an uncertainty that was fitted on our data with a 
linear regression (Fig. 4-b) 



cr D = 0.4525D - 99.32, 



(8) 



which is an upper limit of the uncertainty on the stellar distance. 

The comparison of the distribution of the stellar distance 
showed that we also had to bias the BGM requests according 
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to the luminosity classes, as could be expected given the sam- 
ple selection criteria outlined in Sec. 2. We preferred to use the 
spectroscopic surface gravity since it is the parameter directly 
derived with MATISSE, with a bin size corresponding to a lu- 
minosity class size, i.e. 1.0 dex 2 . We therefore constrained the 
BGM distributions of magnitude J, colour (J - K„) and surface 
gravity to be identical to the observed ones. This indeed makes 
the distance distributions more compatible. The resulting distri- 
butions are presented in Fig. 9. We refer to this simulated sample 
as BGM1 sample hereafter. The agreement is good in all three di- 
rections, but it seems that we observed closer stars than predicted 
by the BGM. We should, however, point out that this assumes a 
correct extinction law, which is difficult to test, particularly in 
the SRcOl direction where differential absorption is noted. The 
observations and BGM1 distance distributions agree within 2<x. 

When building the metallicity distributions for our samples 
and the BGM one (see Fig. 10), we noted that the shape and 
agreement with the BGM1 distributions depend strongly on the 
biasing in luminosity class, hence on the distance distribution 
of the observations. Subsequently, our metallicity distributions 
agree with the BGM1 distributions within 3<x and this is strongly 
correlated with the distance distribution. We conclude from these 
comparisons that it is indeed mandatory to simulate properly the 
selection biases according to the luminosity classes, since they 
have a direct impact on the observed distance. 

3.3. Velocity components 

To build statistically relevant distributions for the velocity com- 
ponents, we applied the same procedure as described above 
for the stellar distance distributions (2 000 Monte-Carlo realisa- 
tions). We propagated the uncertainty on the distance described 
in section 2.3 and we also assumed a typical uncertainty of 4 mas 
year - 1 , and 0.3 km s~ 1 , for the proper motion and radial velocity, 
respectively, in agreement with the values found for our sam- 
ple (see section 2.4). We also performed other requests to the 
Besancon Galatic model towards the different studied directions 
taking the uncertainties we have in our data into account . Then 
we applied the same procedure to build the distribution of kine- 
matical parameters for BGM 1 simulated stars without adjusting 
any bias of the BGM simulations. 

Figure 1 1 illustrates the comparison of the velocity compo- 
nents for our observations with the BGM1 simulated data. For 
each distribution in Fig. 11, we adjusted a Gaussian function 
and reported the mean value and standard deviation in Table 5. 
The shape of BGM1 and the observed velocity components dis- 
tributions agree within 3<x. Few cases of discrepancy larger than 
3<r can, however, be noted. One source of these discrepancies is 
the distance scale that does not agree within 3cr and which has a 
direct impact on the W velocity component. This effect is indeed 
slightly alleviated in the SRcOl direction where the agreement on 
the stellar distances is better for the closest stars, which represent 
the majority of this sample. The other noticeable difference is for 
the V component in the SRcOl direction. It could be due to little 
statistics numbers in this direction. Velocity dispersions are in 
good agreement for the three directions. 

As a conclusion, the agreement (within 3cr) between our de- 
rived kinematic parameters and those simulated with the BGM1 
validates our approach and the results we achieved on the ob- 
served sample. The small differences we pointed out might 
be due to the presence of stellar populations that could differ 



from those assumed when building the observed and simulated 
(BGM1) samples but also to the stellar extinction laws we used 
because they were rather poorly constrained. 

4. Identifying the stellar populations in the 
CoRoT fields 

Depending on the observed Galactic direction, any line of sight 
might contain a mixing of the different stellar populations (thin 
disc, thick disc, halo, and bulge). To properly interpret the ob- 
served data, we need to identify and analyse these different 
stellar populations separately. Figure 12-a) shows the height 
to the Galactic plane of symmetry (Z gc ) as a function of the 

Galactocentric radius (Rcc = yj^cc + ^cc-*- Since our observa- 
tions are very close to the Galactic plane, we expect our sample 
to be mainly composed of thin disc stars. However, the LRcOl di- 
rection is very likely to be contaminated by other Galactic popu- 
lations, since it goes down to 6 kpc below the Galactic plane (see 
Fig. 12-a). We also hypothesised in the previous section that the 
differences between the BGM1 and the observed velocity com- 
ponents distributions could be due to a different mix of stellar 
populations. To confirm this, we took advantage of the BGM1 
simulations described in the previous section. We used the age 
flag given by the BGM1 to differentiate thin disc, thick disc, and 
halo stars (see Tab. 6). The results are summarised in Table 7 
(the first three lines with the age flag). According to the BGM1, 
two of the CoRoT fields are mainly composed of thin-disc stars 
(~95%) but the LRcOl direction indeed contains a non negligi- 
ble amount (17%) of thick-disc stars. None of the fields should 
contain halo stars following BGM1. 



Table 6. Age criteria associated with the different stellar popu- 
lations in the BGM model 



Thin disc 


to 10 Gyr 


Thick disc 


> 10 Gyr 


Halo 





One independent way to identify the stellar populations is to 
apply the procedure described by Bensby et al. (2005) to BGM1 
simulations. This procedure combines the three Galactic veloc- 
ity components to build a probability for a given star to belong 
to a given population. We checked that the Galactic rotation re- 
mains close to the solar one for the mean distances covered by 
our sample. We used only the three basic Galactic components 
(thin disc, thick disc, and halo) and we assessed those belonging 
to one population if the probability (defined in Eq. 9) is ten times 
higher than the two others, which is a very strict criterion. The 
results shown in Table 7 (three middle lines), are consistent with 
the age flag given by BGM, although fewer thick disc stars in the 
LRcOl direction are identified. 

The adopted kinematic criteria are (Bensby et al. 2005): 



P = X -k ■ exp 



U 2 

'24 



(V] 



LSR 



V "I 2 
* asym/ 



2a% 



W 2 ^ 

LSR 



where P is either thinD, thickD, or H, and 
1 



(9) 



(10) 



2 We note that different gravity bins have been tested (0.75, 1.25, 
1.5 dex) with no significant changes in the results that are presented. 



(27r) 3 / 2 CF\j 0" V a"w 

is the normalisation factor; cry, cry, cr w are the characteristics 
dispersion velocities; V asym is the asymmetric drift; and X is the 
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Table 7. Galactic population repartition in the BGM1 biased 
simulations. 





Thin disc Thick disc 


Halo 


LKaUl riLjMt age 


96% 


4% 


U7C 


LRcOl BGM1 age 


82% 


17% 


1% 


SRcOl BGM1 age 


95% 


5% 


0% 


LRaOl BGM1 kine. 


87.7% 


2.6% 


0.0% 


LRcOl BGM1 kine. 


76.0% 


6.2% 


0.5% 


SRcOl BGM1 kine. 


86.1% 


3.6% 


0.0% 


LRaOl OBS. kine. 


58% 


18% 


0.3% 


LRcOl OBS. kine. 


48% 


20% 


5% 


SRcOl OBS. kine. 


82% 


2.4% 


0.0% 



Notes, "age" or "kine" correspond to the method used to separate the 
stellar populations: "age" is related to the age flag returned by the 
BGM1 and "kine" denotes classifications based on the Bensby et al. 
(2005) method using velocity components. For the kinematical results, 
the sum of the three component does not reach 100% because it is a 
probabilistic way to disentangle the stellar populations, and some stars 
have intermediate kinematical parameters. 

fraction of observed stars in the solar neighbourhood belonging 
to each population. 

From the BGM age information, we search for misclassified 
stars using kinematical criteria. Indeed our identification of the 
stellar populations depends on kinematics criteria alone, and not 
on chemical information. Besides this, the procedure is based 
on solar neighbourhood data, so by applying it to our stars, we 
assume that these kinematical properties are also valid far from 
solar neighbourhood (Bensby et al. 2011b). Bensby & Feltzing 
(2010) show that these kinematical criteria could introduce sig- 
nificant mixing in the identified stellar populations, but it has the 
advantage of not imposing an a priori on their chemical proper- 
ties. 

Applying these kinematics criteria to the BGM1 simulations, 
we found that thin disc stars are not contaminated by any of 
the two other populations. However, about 20% of the simulated 
stars classified as thick disc members in the LRaOl direction are 
actually thin disc stars. This contamination of the thick disc pop- 
ulation by thin disc stars is negligible (~1%) in the SRcOl direc- 
tion and goes up to ~50% in the LRcOl direction. Finally, the 
halo population actually comprises ~10% of thick disc stars and 
~20% of thin disc stars in the LRcOl direction. The halo is very 
poorly represented in all of the simulated samples, showing that 
this population can be ignored altogether in the interpretation of 
our observed samples. For that reason, we decided in what fol- 
lows to only characterise the chemo-dynamical properties of the 
thin disc, since obtaining a clean sample of the other components 
is questionable and could lead to misleading results. In addition, 

Table 8. Characteristic velocity dispersion (<t\j, <x v , and crw) of 
the thin, thick disc, and halo for the equation (9). V asym is the 
asymmetric drift and X the fraction we chose. 







Cy cr w 


^asym 


X 






(kms- 1 ) 




% 


thin disc (thinD) 


35 


20 16 


-15 


94 


thick disc (thickD) 


67 


38 35 


-46 


6 


halo (H) 


160 


90 90 


-220 


0.15 



we stress that, whereas the BGM assumes a variation in the kine- 
matics with the galactocentric distance, no particular adjustment 
of the velocity ellipsoid is needed for the spanned distances to re- 
trieve the properties of each simulated population. Nevertheless, 
imposing the same kinematics as in the solar neighbourhood for 
the thick disc and the halo at such distances might not be realis- 
tic, so would at least bias our results kinematically. This is thus 
another reason for not analysing these two galactic structures. 

Finally, to measure the properties of the thin disc, we used 
a simulation of 500 Monte-Carlo realisations, for which the ob- 
served stellar velocities were drawn from a normal distribution 
around their measured value and with a standard deviation equal 
to the assigned uncertainty on the measurement. For each of 
the realisations, the kinematical criteria have been applied and 
a membership for each star has been associated, as described 
previously (Fig. 12-b to e, show the results for one particular 
realisation: the one with the exact measured values). The mean 
kinematics, gradients and metallicities were then estimated for 
each realisation, by fitting a Gaussian to the derived histograms. 
The values presented in Table 7 (bottom three lines) represent 
the mean proportion of each population after the Monte-Carlo 
simulation. As expected, the three directions are mainly com- 
posed of thin disc stars. The SRcOl direction is the most con- 
sistent one with the BGM1 because it is almost entirely com- 
posed of thin disc stars. As a result, the differences observed 
in the stellar distance and the velocity component distributions 
for the SRcOl are probably due to the assumed extinction laws, 
and the small statistics numbers. The LRaOl direction contains 
more thick disc stars than predicted by the BGM1, which could 
be one of the sources of the differences between the BGM1 and 
our observations. These two directions (SRcOl and LRaOl) are 
therefore both mainly composed of thin disc stars with a very 
negligible amount of thick disc and no halo stars. Finally, the 
LRcOl field, however, contains a weaker thin disc contribution. 
It is the most mixed one and may be slightly richer in terms 
of thick disc stars than expected by the BGM1 simulations, al- 
though we cannot exclude that this excess thick disc population 
is made up of misclassified thin disc stars. Indeed, as pointed out 
by the analysis of the BGM1 simulated stars, although little er- 
ror on the classification is expected for stars attributed to the thin 
disc population, the misclassification is expected to be higher for 
the thick disc stars, with a significant fraction of them being thin 
disc stars. It is also probable that there are very few thick disc 
stars misclassified as thin disc ones since our targeted galactic 
directions contain fewer thick disc members. 

All these remarks show that it would be safer to disentangle 
these stellar populations before any study devoted to the thin disc 
in these fields. These results are expected from the respective 
mean Z gc and the distance range of these fields (see Fig 12-a) and 
could explain the small differences between the space velocity 
components illustrated in Fig. 9, 1 1 and Table 5. 

5. Properties of the thin disc 

According to the strict kinematics criteria and the Monte- 
Carlo simulations described in the previous section (bottom of 
Table 7), there are on average 240 thin disc stars out of 404 stars 
observed in the LRaOl direction (-59%), 140 out of 286 (-49%) 
in the LRcOl direction, and 52 out of 64 (~8 1 %) in the SRcOl di- 
rection. In this section, we discuss the kinematical and chemical 
properties of these stars identified as thin disc population. When 
analysing the kinematics of the separated stellar populations in 
the three fields, one has to keep in mind that we used this infor- 
mation to select the stars belonging to the thin disc. We recall 
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Table 9. 

disc 



Kinematical and chemical characteristics of the thin 











< N t > 


240 + 6 


140 ±6 


52 ±2 


<U> 


-24.4 ± 1.4 


4.6 ±2.3 


11.9 ±3.3 


<v> 


-18.6 ± 1.4 


-13.0 ± 2.3 


-7.4 ± 2.6 


< w > 


-0.9 ± 1.5 


-6.2 ±2.1 


-7.8 ±5.1 


CTu 


27.8 ± 1.9 


30.3 ±3.1 


29.1 ±4.5 


<j v 


13.0 ±3.2 


21.7 ±4.3 


13.1 ±4.5 


<T W 


17.2 ±6.8 


15.8 ±6.9 


12.2 ±4.8 


<[M/H]> 


-0.05 ± 0.01 


0.07 ± 0.02 


0.07 ± 0.02 


°~ [M/H] 


0.23 ± 0.01 


0.28 ± 0.02 


0.23 ± 0.04 


<[tt/Fe]> 


0.06 ± 0.01 


0.04 ± 0.01 


0.05 ± 0.01 


C [ff/Fe] 


0.05 ± 0.01 


0.07 ± 0.01 


0.05 ± 0.01 



Notes. The solar motion have not been corrected here. The dispersion 
values, for each Monte-Carlo realisation, have all been corrected by the 
observational errors as in Jones & Walker (1988) 
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LRc01 
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Fig. 13. Distribution of the overall metallicity for the thin disc 
sample in the three pointing directions. 



that, based on the above discussion, it is probable that this thin 
disc sample is relatively pure, although biased at some level. 

5.1. Kinematics 

By selecting the thin disc stars, we ended up with the closest 
stars. The most distant thin disc stars are found in the LRcOl 
at Zqq — -350 pc, which is compatible with what is generally 
admitted for the scale of height of the thin disc (see Fig. 12-a 
and Veltz et al. 2008). We used the same method as described 
in Sect. 4 to build the distributions for each Monte-Carlo reali- 
sation and to estimate the means and standard deviations of the 
velocity components in each direction. The dispersions reported 
in Table 9 are corrected from the observational errors as in Jones 
& Walker (1988), for each realisation. As expected, the values 
in Table 9 are in generally good agreement with the criteria we 
used for disentangling the stellar populations. In addition, the ve- 
locity dispersions appear to correspond to the oldest part of the 
thin disc simulated in the BGM, though for the V-component we 
probably underestimate the dispersion due to an overestimation 
of the individual error on the V measurements. Finally, we anal- 
ysed the evolution of the mean V-velocity for the entire sample, 
as a function of the galactocentric radius and found a gradient of 
-5.3 ± 1.3 km s" 1 kpc _1 . 

5.2. Metallicity 

Our entire thin disc sample contains stars with metallicity rang- 
ing from -0.88 ± 0.17 to 0.55 ± 0.20 dex. The metallicity dis- 
tributions peak around 0.07 dex for the LRcOl and SRcOl di- 
rections close to the Galactic centre and around -0.05 dex for 
the LRaOl direction, with a standard deviation of ~0.25 dex (see 
Fig. 13). The mean values reported in Table 9 agree, within lcr, 
with the mean values used for the thin disc simulation in the 
BGM (Robin et al. 2003) as expected according to the observa- 
tional uncertainties. The LRaOl field exhibits a mean metallicity 
completely in agreement with the values for the thin disc pub- 
lished by Fuhrmann (2004) and for the SDSS sample presented 
by Allende Prieto (2010). However, the fields in the centre di- 
rection are slightly more metal rich. 



5.3. Radial mixing and metallicity gradient 

In Fig. 13, stars located towards the Galactic centre appear to 
be more metal rich than those found in the anti-centre direction. 
In Paper I, we separated the giant from the dwarf stars using 
the measured log g and noted that this effect is clearer for giant 
stars, probably because the giant stars may be observed further 
than dwarfs. Here we use the same criterion to separate the giants 
(log g< 3.5 dex) from dwarfs (log g> 3.5 dex). These metallicity 
distributions can be different because of the radial and vertical 
metallicity gradients, radial mixing processes, and local inho- 
mogeneities (Pedicelli et al. 2009; Schonrich & Binney 2009; 
Haywood 2008). We propose here to search for radial mixing 
signatures and to quantify the radial metallicity gradient. 

Navarro et al. (201 1) find no correlation between V velocity 
and metallicity, thereby constraining the degree of radial mix- 
ing (Schonrich & Binney 2009; Haywood 2008). We searched 
for such a correlation in our thin disc sample. For each of our 
Monte-Carlo realisations, we have measured the correlation 3 be- 
tween the metallicity and the V-velocity, by selecting stars in 
narrow ranges of galactocentric radii. Indeed, evidence of radial 
migration is only visible for localised samples. When larger dis- 
tances are involved, correlation between metallicity and velocity 
appear naturally, just because the metal-richer stars in the in- 
ner galaxy rotate slower, due to the higher vertical velocity dis- 
persion in these regions. We considered both giants and dwarf 
stars in the radial ranges 7 < Rcc < 8 kpc, 8 < Rcc < 9 kpc 
and 9 < Rcc < 10 kpc, and find the following correlations 
R = -0.028 + 0.077, -0.094 ± 0.052, -0.192 ± 0.059, respec- 
tively. Whereas for the most inner part the correlation is not sig- 
nificant considering the uncertainties, this is not the case for the 
two other distance bins. The significance of the value at ~ 2<x 
is characteristic of radial migration processes in the thin disc. 
Indeed, as a sanity check we searched for such a correlation in 
the BGM1 simulation, since we recall that the BGM does not in- 
corporate any radial migration. No such correlation was found, 
as expected, which is consistent with radial-mixing processes in 
our observations. 

The thin disc metallicity gradient is also an important in- 
put for stellar populations models to understand galactic for- 
mation and evolution (Schonrich & Binney 2009). The value 
of the metallicity gradient is a matter of debate, ranging from 
-0.04 to -0.1 dex kpc 1 (Haywood 2008). We explore different 

3 Spearman's rank correlation value 
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Galactic regions, within 1.4 kpc and 600 pc, with the giants and 
the dwarfs, respectively. We calculated the weighted 4 average of 
the Galactocentric radius and the metallicity for the LRaOl field, 
on the one hand, and the LRcOl and SRcOl fields on the other. 
We found for the dwarfs a gradient of -0.053 ± 0.01 1 dex kpc -1 , 
and for the giants, -0.097 ± 0.015 dex kpc -1 . These results are 
illustrated in Figs. 14-a) and b) for the dwarfs and for the gi- 
ants, respectively. On the one hand, the dwarf stars are known 
to have better determined atmospheric parameters, hence better 
distance and metallicity. However, the distances covered by the 
giants range further than for the dwarfs, which allow us to ex- 
plore a larger difference in metallicity and a larger differential 
of distance. On the other hand, for giant stars, the uncertainties 
on the log g might result in underestimating the distance of the 
stars from the Sun, thus overestimating the metallicity gradient. 
We note, however, that when limiting it to the same distance 
range, the measurement of the gradient for the dwarfs and the 
giants marginally agree, though the same trend is apparent. For 
the same radial range between 7.1 kpc and 9.3 kpc, we mea- 
sure gradients of -0.041 ± 0.013 dex kpc -1 for the dwarfs and 
-0.057 ± 0.019 dex kpc -1 for the giants. 

On the other hand, we point out that, in order to test possi- 
ble effects due to some bias in our procedure, we measured the 
metallicity gradient in the BGM1 sample and showed that the re- 
sult with and without uncertainties on kinematical and chemical 
parameters agree within error bars. 

We also tried to estimate the effect of our sample selection 
(magnitude limited sample) and individual distance uncertain- 
ties when deriving average quantities in discrete distance bins 
(Malmquist-type biases) such as the gradients. For that pur- 
pose, we derived the gradients from Besancon simulated stars 
by adopting the whole sample and by randomly rejecting some 
faint stars (i.e. by mimicking the spectroscopic selection). It has 
been found that the results are only weakly affected by this ef- 
fect. Finally, we compared the metallicity gradients estimated for 
the BGM and BGM1 samples (direct test of the Malmquist bias). 
These gradients differ by ~20% only if we consider all the stars 
found in BGM and BGM1. Furthermore, the difference between 
the BGM and BGM1 gradients is around ~35% if we only take 
the giant stars into account. The BGM1 gradient is always the 
steepest one in these different tests. This is an indication that the 
Malmquist bias could affect our estimated metallicity gradients, 
and these gradients could be slightly overestimated. However, 
their values are still reasonably well determined, with potential 
biasses that do not exceed the quoted uncertainties (see also their 
comparison with other literature determinations below). 

The measured metallicity gradient using dwarf stars agree, 
within the error bars, with the most recent studies using Cepheid 
stars, with Galactocentric distances ranging from ~5 to ~17 kpc 
(-0.051 ± 0.004 dex kpc -1 , Pedicelli et al. 2009), and slightly 
higher when we only consider the giants. Moreover, these au- 
thors show that the gradient is much steeper for the inner Galaxy 
(-0.13 ± 0.015 dex kpc -1 ). We were not able to separate the 
inner and outer Galaxy in our sample since the scatter of each 
subsample is too high. Our metallicity gradient value is also 
slightly higher but nevertheless compatible within 1 cr with stud- 
ies based on Hn regions. The latest work by Rudolph et al. 
(2006) explores regions with Rqc ranging from 10 kpc to 15 kpc 
and presents metallicity gradients 5 ranging from -0.046 ± 0.009 
to -0.071 + 0.010 dex kpc -1 using optical data. These authors 
showed that far-infrared data results in different gradient val- 
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Fig. 14. a) Overall metallicity as a function of Galactocentric ra- 
dius for dwarf stars in the thin disc, b) Overall metallicity as a 
function of Galactocentric radius for giant stars in the thin disc, 
c) ^-enhancement as a function of Galactocentric radius. We 
plotted the gradient found for giant (dash-dotted line) and dwarf 
(dashed line) stars. Blue X represent the stars from LRcOl, green 
circles represent stars from SRcOl, and red + represent stars from 
LRaOl. 

ues ranging from -0.041 + 0.009 to -0.085 ± 0.010 dex kpc -1 . 
They finally emphasise that these results are extremely sensitive 
to extinction, which is one of the problems we had to face in 
our own study. Finally, observations of open clusters provide a 
fair estimate of metallicity and Galactocentric radius for older 
stars. Such studies find steeper gradients, in very good agree- 
ment with our estimate for the giant stars (Andreuzzi et al. 201 1 ; 
Magrini et al. 2010; Carraro et al. 2007; Friel et al. 2002). We 
also note that our radial metallicity gradient value for the gi- 
ant stars, probing the furthest regions of the Galactic plane, may 
be compatible with the radial gradient used in the BGM. It is 
also fairly compatible with what is found using other methods. 
Similar studies are also devoted to other nearby galaxies such as 
M81, for which Stanghellini et al. (2010) find an oxygen gradi- 
ent of -0.055 ± 0.02 dex kpc -1 , close to the value we derived for 
the Milky Way dwarfs. 

5.4. a-enhancement 

The measurement of the [c/Fe] ratio for different stellar popula- 
tions is important for our understanding of chemical history of 
the Galaxy. Table 9 reports the mean and cr values for the three 
fields. The mean [c/Fe] is identical in the three Galactic direc- 
tions, and compatible with thin disc canonical value. We used 
the same method as described above to measure the radial gradi- 
ent for [a/Fe]. We found no variation in [c/Fe] as a function of 
the Galactocentric radius (see Fig. 14-b). 

6. Stars with unusually high [a/Fe] 

Among our sample of 754 stars, we found 110 stars with 
[ff/Fe] values substantially higher, i.e. beyond the error bars 
(o"[a/Fe] -0.1 dex), than the standard description based on galac- 
tic stellar properties (see Paper I, and solid line in Fig. 12-f 
where they can be easily identified). These stars are found above 
+ lcr and none are found below -lcr, which makes these de- 
tections statistically significant. The majority of these stars are 
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main sequence stars with mean metallicities in the range -0.5 
to 0.5 dex. There is, however, a greater proportion of giant stars 
among this group than in the entire sample. It is also instruc- 
tive that ~70% of these stars are located in the LRcOl direc- 
tion. These stars show no peculiarity in terms of velocity com- 
ponents, and they are found at any distance towards the centre 
direction, and within the closer kiloparsec towards the LRaOl 
direction. They mainly belong to the thin disc (56%). We per- 
formed some additional tests to further verify the MATISSE pa- 
rameters and particularly the [a/Fe] values. We first checked for 
these stars that the T e ff derived by MATISSE are consistent with 
their observed colours in the infrared (using 2MASS photom- 
etry). We then selected those among these stars with the best 
spectra (SNR> 20, i.e. 28 stars). We performed a visual check 
for these stars, and we also searched for the synthetic spectrum 
in the grid that is the most similar to the observed one by min- 
imising the x 2 between observed and synthetic spectra. This al- 
lowed us to confirm the MATISSE parameters for eleven of these 
stars. Another ten of these spectra present a slightly better fit 
with atmospheric parameters that re different from MATISSE 
ones, but the difference never exceeds a step of the grid and are 
still pointing to a-enhanced abundances. The seven remaining 
spectra are not adjusted with the parameters of MATISSE, nor 
with those of the fit based on the minimum^ 2 . The examina- 
tion of the cross-correlation function for these seven spectra has 
revealed that these stars might be spectroscopic binaries not de- 
tected by Loeillet et al. (2008). Therefore from these 28 rather 
good SNR spectra, we can confirm the peculiar [a/Fe] for 21 
stars. Finally, we extended the x 2 test to all of these ^-peculiar 
stars: we searched for the grid point minimising the x 1 with the 
observed spectrum. We found that 90% of these spectra have x 2 
parameters a-enhanced. 

All these checks therefore confirm the non-standard 
[a/Fe] value of this subsample of stars. We therefore might 
have detected metal-rich stars with unexpected peculiar high 
[a/Fe] values that contradict the current scenarios of evolution 
and chemistry of the galactic disc. For example, it has been pro- 
posed that bulge stars could be a-enhanced because they are 
supposed to form quickly (Mc William & Rich 1994; Matteucci 
& Brocato 1990), although more recent studies seem to con- 
tradict this point (see, for instance, Bensby et al. (2011a)). Our 
a-peculiar stars might find an origin in the Galactic bulge and 
have migrated to the solar neighbourhood via radial mixing. 
These stars would deserve further studies to estimate their chem- 
ical abundances in detail and to disentangle their origin. 6 

7. Properties of the planet-hosting stars in the 
targeted directions 

In Paper I, we showed that the de-biased metallicity distributions 
combined with planet-occurrence probability laws, giant-planet 
period distribution, and geometric probability of transit, pro- 
vided a number of detections agree with what was detected by 
CoRoT in these fields. CoRoT detected seven planetary systems 
in the Galactic directions targetted in this study. We retrieved 
the atmospheric parameters from the literature (see Table 10), 
except for CoRoT-lb for which we had already derived param- 
eters in Paper I. These atmospheric parameters agree, within 
the error bars, with the values published in Leger et al. (2009) 
and Bruntt et al. (2010). From these parameters, we used the 

6 A few months after the submission of this article, a discovery of 
similar a-enhanced metal-rich stars has been reported by Adibekyan 
et al. (201 1), thus confirming our detections. 



same methodology as for our GIRAFFE/FLAMES sample to de- 
rive the kinematical information. The distance derived here is 
generally in good agreement with published values for these 
planet-hosting stars. We represented these planet-hosting stars 
in Fig. 12-a). We also applied the same kinematical criteria to 
identify the stellar population these stars might belong to. We 
find that the planet-hosting stars mainly belong to the thin disc. 
However, the kinematics of CoRoT-3b and 8b did not allow us 
to classify them if we adopt the same strict criterion as for our 
GIRAFFE/FLAMES sample. Indeed, the probability of belong- 
ing to the thin disc is only 2.5 times greater than the probability 
of belonging to the thick disc for CoRo T-3b and only 0.5 times 
greater for CoRo T-8b. It is particularly interesting to note that 
these two stars are located in the LRcOl field, which is the field 
that is the most mixed up in terms of populations and ages ac- 
cording to BGM. However, the Zee for these stars is compatible 
with thin disc stars, as shown in Fig. 12-a). We also find a clear 
correlation (R=-0.71) between the V velocity component and 
the metallicity for these planet-hosting stars. Although the statis- 
tics is low, this tends to favour a radial mixing origin of planet 
metallicity correlation as proposed by Haywood (2009). 

8. Conclusion 

We computed stellar distances for the 754 CoRoT stars with 
atmospheric parameters derived using the MATISSE algorithm 
(Gazzano et al. 2010). This allowed us to map the Galactic kine- 
matics and chemistry in three Galactic directions observed by 
CoRoT. All the results are available electronically through the 
Exo-Dat database (Deleuil et al. 2009). 

Using the kinematical criteria described by Bensby et al. 
(2005), we identified thin disc and thick disc stars in our ob- 
served sample. The procedure was first validated on stars sim- 
ulated with the Besancon Galactic model which allowed us to 
investigate the possible limitation of the classification method. 
We found that the proportion of thin disc stars misclassified as 
thick disc ones should be significant especially in the LRcOl . Our 
results for the observed sample show that our kinematical prop- 
erties are in good agreement with predictions of the Besancon 
model, which is, however, very sensitive to the stellar extinction 
law assumed. The LRaOl and LRcOl directions contain fewer 
thin disc stars and more stars suspected of belonging to the thick 
disc population than predicted. This could explain the small dif- 
ferences observed between space-velocity component distribu- 
tions. 

The adopted selection criteria also allowed us to build a 
clean sample of thin disc stars to study the properties of this 
stellar population. Combining the velocity dispersions and the 
chemistry, we found a correlation between the V-velocity com- 
ponent and [M/H], for two of the three considered radial bins, 
which could be a clue of radial mixing (Navarro et al. 2011; 
Schonrich & Binney 2009; Haywood 2008). With dwarf stars, 
we also measured a thin disc radial metallicity gradient of 
-0.053 ± 0.011 dex kpc~', which is consistent with the most 
recent published values (Maciel & Costa 2010). 

Our analysis also shows the presence of stars in our sample 
with unexpected high [a/Fe] values at rather high metallicities. 
This stellar population might have been formed in the Galactic 
bulge and migrated up to solar neigbourhood (McWilliam & 
Rich 1994; Matteucci & Brocato 1990). Although our spectro- 
scopic data are not appropriate for carrying out the detailed anal- 
ysis that such stars would require to better understand their ori- 
gin, the existence of this puzzling population has nevertheless 
been recently confirmed by Adibekyan et al. (201 1). 
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Applying the same methodology to the planet-hosting stars 
discovered by CoRoT in the targeted directions, we showed that 
they mainly belong to the thin disc stellar population. This prob- 
ably means that, at least in these three fields, CoRoT has detected 
planets only around thin disc stars. We also found a correlation 
(R=-0.71) between the V velocity component and the metal- 
licity for these planet hosting stars, suggesting a radial mixing 
history for these stars. 

This work is one of the first studies of the stellar pop- 
ulations in the Galactic plane not limited to the strict so- 
lar neighbourhood, and based on good statistics although still 
limited. This demonstrated the potential of multi-fibre instru- 
ment like FLAMES/GIRAFFE, combined with automatic anal- 
ysis tools like MATISSE for Galactic physics analyses. It would 
be interesting to complete this study in other Galactic direc- 
tions, observed or not by CoRoT, and to combine these results 
with the richness of the information that can be derived from 
CoRoT light-curves analysis. For instance, one could combine 
the atmospheric parameters with the light curve parameters, i.e. 
degree and scale of variability, to search for links between the 
stellar parameters and the photometric variations of the star. 
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Table 2. Un-reddened colours, absorption, absolute magnitudes and stellar distance of the 754 CoRoT stars of the present study 
matching our quality criteria. 



CoRoT JD (J - K s ) 


A(7 - K s ) 


(V-J) A(V-J) 


Aj 


AAj 


My 




D(pc) 


AD (pc) 


211652185 0.624 
211652936 0.367 
211666039 0.296 


0.233 
0.099 
0.096 


1.738 0.231 
1.183 0.164 
1.028 0.145 


0.812 
0.143 
0.315 


0.412 
0.174 
0.183 


1.550 
4.640 
3.630 


0.610 
0.460 
0.690 


905 
357 
511 


321 
85 
171 


Table 3. Kinematics results for the 754 stars fulfilling our quality criteria. 


C0R0T.ID m n b 
(mas cent ) 


X AX 


Y AY Z 
(kpc) 


AZ 


V 


AU 


V 

(kms~ 


AV 

- 1 ) 


W AW 


211652185 -0.27413 -1.24966 
211652936 13.30323 4.81142 
211666039 -17.94455 -7.38882 


-7.781 0.255 
-8.217 0.068 
-8.088 0.138 


0.551 0.196 -0.017 
0.218 0.052 -0.007 
0.302 0.101 -0.017 


0.006 
0.002 
0.006 


-3.3 
-34.0 
46.8 


12.3 
5.1 
10.0 


-4.0 
2.3 
-19.6 


16.0 
6.7 

13.9 - 


-5.3 20.3 
8.6 6.7 
18.8 10.8 




7500 7000 6500 6000 5500 5000 4500 4000 7500 7000 6500 6000 5500 5000 4500 4000 

Teff [K] Teff [K] 

Fig. 6. Hertzsprung Russel diagram in the T e ff-log g plane for 2 000 random entries of the BGM request in the LRcOl direction (left) 
and the resulting MATISSE atmospheric parameters derived from the corresponding synthetic spectra at SNR=10. 



Table 5. Kinematical results for the observations and the BGM1 biased simulations. The unit is km s 





<U> 


cr v 


< V > 


cr v 


< W > 


0"w 


LRaOl 


-33.6+1.5 


40.3+1.2 


-20.0+1.5 


34.7±1.3 


1.3+1.6 


35.1+1.4 


LRaOl bgm 1 


-27.8±0.2 


35.0±0.1 


-16.6±0.2 


33.8±0.2 


-5.5±0.2 


31.4+0.2 


LRcOl 


16.7+2.1 


50.9±1.7 


-17.5±2.0 


45.2±1.6 


-13.5+2.1 


45.8+1.8 


LRcOIbqmi 


21.1+0.1 


48.4±0.1 


-17.3+0.1 


38.8±0.1 


-5.3+0.1 


36.9±0.1 


SRcOl 


11.9+2.6 


42.2±1.6 


-6.8±1.9 


25.6+1.5 


-7.8±1.9 


22.4±1.5 


SRcOIbqmi 


10.1+0.1 


40.1±0.0 


-12.9±0.0 


27.5±0.0 


-7.0±0.0 


24.5±0.0 



Table 10. Atmospheric parameters for the planet-hosting stars detected by CoRoT in our targetted directions 



CoRoT field 


CoRoT ID 


Planet Name 


Teff (K) 


°T eff 


logg 


°"i°g g 


[M/H] 


C [M/H] 




O" [a/Fe] 


Reference 


LRcOl 


101206560 


CoRoT-2b 


5625.0 


163.5 


4.3 


0.335 


0.0 


0.182 


0.0 


0.217 


(1) 


LRcOl 


101368192 


CoRoT-3b 


6740.0 


178.7 


4.22 


0.278 


-0.02 


0.163 


-0.05 


0.115 


(2) 


IRaOl 


102912369 


CoRoT-4b 


6190.0 


126.2 


4.41 


0.274 


0.05 


0.167 


0.0 


0.217 


(3) 


LRaOl 


102708694 


CoRoT-7b 


5319.0 


120.4 


4.76 


0.186 


0.21 


0.17 


0.05 


0.071 


(4) 


LRcOl 


101086161 


CoRoT-8b 


5080.0 


143.4 


4.58 


0.201 


0.3 


0.196 


-0.04 


0.099 


(5) 


LRcOl 


100725706 


CoRoT-lOb 


5075.0 


140.7 


4.65 


0.209 


0.26 


0.196 


0.06 


0.139 


(6) 


LRaOl 


102671819 


CoRoT-12b 


5675.0 


136.8 


4.52 


0.281 


0.16 


0.182 


-0.06 


0.138 


(7) 



References. (1) Alonso et al. (2008); Bouchy et al. (2008) ; (2) Deleuil et al. (2008) ; (3) Moutou et al. (2008) ; (4) Gazzano et al. (2010) ; (5) 
Borde et al. (2010) ; (6) Bonomo et al. (2010) ; (7) Gillon et al. (2010) 
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Fig. 7. Evolution of the bias (black crosses) and dispersion (red asterisks) of the MATISSE atmospheric parameters for synthetic 
spectra (at SNR=10) of 2000 stars randomly chosen in the LRcOl BGM request. 
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Fig. 9. Distribution of the stellar distances for our observed sample (red dashed line) and the BGM1 simulation (blue). The error 
bars correspond to 1 <x of the Gaussian function fitting the distribution of the corresponding bin. 
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Fig. 10. Distribution of the metallicity for our observed sample (red dashed line) and the BGM1 simulation (blue). The error bars 
correspond to lcr of the Gaussian function fitting the distribution of the corresponding bin. 
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Fig. 11. Kinematic comparison (three velocity components) of our observed sample with the BGM1 simulations for the three 
CoRoT fields. The observations are represented with a red dashed line and the BGM1 data with a solid blue line. 



Table 11. Absorption, distance, and geocentric coordinates for the CoRoT planet-hosting stars in our fields 



CoRoT ID 


Aj 


<J~Aj 


D 




X 




Y 


(Ty 


Z 


""z 








(pc) 








(kpc) 






CoRoT-2b 


0.1436 


0.2047 


275.387 


84.718 


0.214 


0.066 


0.169 


0.052 


-0.035 


0.011 


CoRoT-3b 


0.2659 


0.1855 


789.243 


231.549 


0.622 


0.182 


0.473 


0.139 


-0.111 


0.033 


CoRoT-4b 


0.1156 


0.1709 


862.254 


224.761 


-0.721 


0.188 


-0.472 


0.123 


-0.014 


0.0040 


CoRoT-7b 


0.0582 


0.1875 


171.124 


28.299 


-0.143 


0.024 


-0.093 


0.015 


-0.0070 


0.0010 


CoRoT-8b 


0.2793 


0.2193 


324.957 


70.942 


0.253 


0.055 


0.199 


0.044 


-0.04 


0.0090 


CoRoT-lOb 


0.3672 


0.2204 


366.143 


77.352 


0.289 


0.061 


0.221 


0.047 


-0.044 


0.0090 


CoRoT-12b 


0.1585 


0.2006 


1144.211 


250.813 


-0.957 


0.21 


-0.625 


0.137 


-0.049 


0.011 
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Fig. 12. a) Height to the Galactic plane of symmetry (Zee) as a function of the Galactic radius (Rcc) f° r the whole sample. Blue 
crosses represent stars from the LRcOl direction, green circles stars from the SRcOl direction, and red plus stars from the LRaOl 
direction. Planet-hosting stars are represented with black thick uncertainty bars. The dotted lines correspond to the typical scale 
height of the thin and thick discs reported by Veltz et al. (2008). b) Toomre diagram representing the different Galactic populations 
in our sample: thin disc (cyan A), thick disc (yellow *) and halo stars, (orange o). Dotted lines represent circles with radius every 
50 kms . c) Bottlinger diagram: U versus V velocities. The symbols are the same as in the Toomre diagram. The thin disc stars 
are focussed at U = kms~' and V — kms -1 . d) Kinematical heat diagram of W velocities. The symbols are the same as in the 
Toomre diagram. The separation between the stellar populations is clear with the kinematically hottest stars being the halo stars and 
the coolest ones corresponding to the thin disc, e) Peculiar velocity (v pec = Vt/ 2 + V 2 + W 2 ) as a function of the overall metallicity. 
The symbols are the same as in (f). We represented in thick black the planet-hosting stars for (a) to (e). f) a-enhancement as a 
function of the metallicity. The solid line corresponds to the standard law we used in our grid, the dashed line corresponds to lcr 
deviation from this law. 
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Table 12. Velocity components for the CoRoT planet-hosting stars in our directions. 



CoKol ID 


U 




V 


(Ty 


w 


cr w 








(km s 








CoRoT-2b 


27.6 


3.31 


2.356 


4.158 


-2.38 


2.467 


CoRoT-3b 


-60.514 


9.77 


-3.056 


14.263 


47.093 


17.883 


CoRoT-4b 


-22.014 


9.278 


-24.873 


14.143 


-16.592 


17.273 


CoRoT-7b 


-22.246 


0.845 


-22.628 


1.366 


7.109 


1.778 


CoRoT-8b 


-9.805 


5.54 


-56.366 


8.562 


-38.62 


11.007 


CoRoT-lOb 


32.593 


5.787 


-20.686 


8.052 


-17.141 


6.993 


CoRoT-12b 


-2.56 


12.293 


-18.959 


18.825 


9.388 


22.324 



